Reaction rates for Neutron Capture Reactions to C-, N- and 
O-isotopes to the neutron rich side of stability 



H. Herndl, R. Hofinger, J. Jank and H. Oberhummer 

Institut fur Kernphysik, Technische Universitdt Wien, Wiedner Hauptstrafie 8-10, A-1040 Wien, 

Austria 

ON. 

J. Gorres, M. Wiescher 
Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556, USA 

< 

F.-K. Thielemann 

Department fur Physik und Astronomie, Universitdt Basel, Klingelbergstr. 82, CH-4056 Basel, 
OO ■ Switzerland 

o 

00 

o . 

Q\ . B.A. Brown 

Department of Physics and Astronomy, Michigan State University, East Lansing, USA 

(February 9, 2008) 

> 

X ■ Abstract 

a: 

The reaction rates of neutron capture reactions on light nuclei are impor- 
tant for reliably simulating nucleosynthesis in a variety of stellar scenarios. 
Neutron capture reaction rates on neutron-rich C-, N-, and O-isotopes are 
calculated in the framework of a hybrid compound and direct capture model. 
The results are tabulated and compared with the results of previous calcula- 
tions as well as with experimental results. 
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I. INTRODUCTION 



Neutron capture processes on neutron-rich C-, N- and O-isotopes play an important 
role in astrophysical scenarios ranging from nucleosynthesis in the stellar helium and carbon 
burning stages to possibly inhomogeneous Big Bang models. To simulate the nucleosynthesis 
of light isotopes between carbon and neon a detailed understanding of the these neutron 
capture reactions is therefore essential. This paper intends to derive a consistent set of stellar 
neutron capture reaction rates for neutron-rich carbon, nitrogen and oxygen isotopes, based 
on the latest experimental and theoretical information. 

In stellar helium core burning in massive Red Giant stars neutrons are abundantly pro- 
duced via the 14 N(a,7) 18 F(/3 + z/) 18 0(a,7) 22 Ne(a,n) reaction sequence 0,pJ. In the early 
helium core evolution 14 N and 18 O are initially depleted, while initiating the neutron pro- 
duction. The neutrons subsequently trigger the weak s-process component which leads to 
the production of intermediate mass nuclei around A=100. Neutrons, however, can also be 
captured on 12 C and 16 O which are abundantly produced in helium burning. 

In the very last phase of helium core burning, the core expands and outer layers with 
high 14 N and 18 abundances cause an increase in neutron production. Neutron capture 
on these isotopes however may act as neutron poison and may also change the light isotope 
abundances. This depends critically on the reaction cross section in the energy range between 
25 and 200 keV. 

He-shell burning in low mass asymptotic giant branch (AGB) stars has been proposed as 
the site for the main component of the s-process p-gj. Neutron production is triggered in 
the He-burning shell by the 13 C(a,n) reaction on 13 C being ingested by convective processes 
during the thermal pulses. It has been suggested that neutron induced processes on the 
initial and the additionally ingested carbon, nitrogen, and oxygen abundances may have 
considerable influence on light isotope nucleosynthesis in the thermal pulse Recent 
stellar model calculations |J indicate that the reaction cross sections for the neutron capture 
processes need to be known in the energy range between 5 keV and 30 keV. 
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In the framework of inhomogeneous Big Bang models (IMs) high neutron flux induces 
primordial nucleosynthesis which bridges the mass 5 and mass 8 gap ||. Subsequent neutron 
capture processes on neutron-rich carbon, nitrogen, and oxygen isotopes may bypass the 
long-lived 14 C and trigger a primordial r-process [p|-|T2|j. The efficiency for the production of 
heavy elements in such a scenario depends sensitively on the respective neutron capture rates 
for these light isotopes. Therefore the neutron capture cross sections have to be determined 
over a wide energy range up to 1.0 Me V. 

Over the last few years considerable effort has been made to determine the neutron 
capture reaction rates for the C-, N-, and O-isotopes experimentally as well as theoretically. 
With the present paper we attempt to summarize the experimentally determined neutron 
capture rates. We present new calculations in the framework of a hybrid compound nucleus 
and direct capture model and compare the results with the experimental data as well as with 
previous calculations using the direct capture model || and the statistical Hauser-Feshbach 



model [p~3| , |14|| . The same models were used to determine in addition neutron capture reaction 
rates on /^"-unstable neutron-rich C-, N-, and O-isotopes. 

In the following Sect. 2 we will present the formalism used in calculating the reaction 
cross sections and reaction rates. We also discuss the experimental and theoretical input 
parameters for the calculations of the different cross sections and reaction rates. In Sect. 3 the 
results of our calculations are compared with previous predictions both with a direct capture 
and a Hauser-Feshbach model. If available we also compare our results with experimental 
data. Finally, in the last section the results are summarized and discussed. 



II. CALCULATION OF THE REACTION RATES 

The cross section for neutron capture processes is dominated by the non-resonant di- 
rect capture (DC) process and by contributions from single resonances which correspond 
to neutron unbound states in the compound nucleus (CN). For calculating the different re- 
action contributions we used a simple hybrid model: the non-resonant contributions were 
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determined by using a direct capture model, the resonant contributions were based on deter- 
mining the resonant Breit-Wigner cross section. In the case of broad resonances interference 
terms have to be taken into account. To determine the neutron capture cross sections on 
the /3-unstable nuclei the necessary input parameter for the calculations (masses, Q-values, 
spin-parity assignments of bound states and resonances, excitation energies, spectroscopic 
factors, density distributions, scattering data) were taken from experimental data. When no 
experimental data were available we used theoretical values, mainly derived from the shell 
model (see below). 

The total reaction rate is given by 

N A < av > tot = N A < av > r +N A < av > nr +N A < av > int , (1) 

where the three terms represent the resonant contribution N A < av > r , the non-resonant 
contribution N A < av > nr and the interferences N A < av >- m f Each contribution will be 
explained in the following sections. 

Another important quantity is the Maxwellian averaged cross section. For a temperature 
kT it is defined by 

< av > kT 2 r°° , . / E 



J(kT) 



J™Ea n>J (E)eiq>(-—yE . (2) 



A. Resonant Reaction Contributions 



The cross section of a single isolated resonance in neutron capture processes is well 



described by the Breit-Wigner formula [15,16 



nh 2 (2 j + i) r n r\ 

2^2(2 Jt + l) {Er -Ef+ (^)' 

where J and j t are the spins of the resonance level and the target nucleus, respectively, E r 
is the resonance energy. The partial widths of the entrance and exit channel are T n and 
T 7 , respectively. The total width r to t is the sum over the partial widths of all channels. 



The neutron partial width T n can be expressed in terms of the single-particle spectroscopic 
factor S and the single-particle width T s p of the resonance state pOlJTl 



r n = c 2 Sxr s . p . , (4) 

where C is the isospin Clebsch-Gordan coefficient. The single-particle width r sp can be 
calculated from the scattering phase shifts of a scattering potential with the potential depth 
being determined by matching the resonance energy. 

The gamma partial widths T 7 are calculated from the electromagnetic reduced transition 
probabilities B(Jj — ► J/;L) which carry the nuclear structure information of the resonance 



states and the final bound states |L7|]. The reduced transition rates were computed within 
the framework of the shell model. 

Most of the transitions in this work are Ml or E2 transitions. For these the relations are 

r E2 [eV] = 8.13 ■ 10" 7 E^[MeV]B(E2)[e 2 fm 4 ] (5) 

and 

r M i[eV] = 1.16 • Kr 2 E3[MeV]B(Ml)[/4] . (6) 
The resonant reaction rate for an isolated narrow resonance can be expressed in terms of 



the resonance strength 077 (in units eV) pO|Jl8 



N A < ov > r = 1.54 • 10 V" 3/2 T 9 3/2 J2( u l)i ex P (-H-605£i/T 9 ) cm 3 mole _1 s _1 , (7) 

% 

where E^ is in MeV and T 9 is the temperature in 10 9 K. The resonance strength wy for a 
resonance is given by 

2i/ -|- 1 r n r«, . , 

The resonance strength has to be determined experimentally by low energy neutron capture 
measurements or has to be derived from the calculated partial widths. 



5 



B. Non-resonant Reaction Contributions 



The non-resonant part of the neutron capture cross section has been calculated using the 
DC model described in P^,BILB2|. The total cross section a nT is determined by the direct 



capture transitions af c to all bound states with the single particle spectroscopic factors 
C 2 Si in the final nucleus 

^ = E(^k DC • (9) 

i 

The DC cross sections af c are determined by the overlap of the scattering wave function in 
the entrance channel, the bound-state wave function in the exit channel and the multipole 
transition-operator. 

In the stellar energy range considered here the non-resonant reaction cross sections are 
predominantly determined by s-wave and p-wave contributions. If the Q- value of the neutron 
capture reaction is clearly higher than the neutron energy the cross section for s-wave neutron 
capture follows the 1/v law. Then the reaction rate is constant over the entire temperature 



range 23]. The s-wave contribution to the reaction rate can then be directly determined 



from the thermal cross section o"th, 

N A < av > s = N A x a th v th cm 3 mole -1 s -1 . (10) 

The cross section for p-wave contributions is approximately proportional to the relative 
velocity v, the reaction rate is therefore proportional to the temperature p| j23|j2^| and can 
be expressed by 

N A <av > p = L08 " 1Q8 .^ffl-T 9 cm 3 mole" 1 s" 1 . (11) 

For low Q- values the simple 1/v- and f-law do not apply anymore. A significant deviation 
can be observed if the neutron energy is in the order of the Q- value. In this case the energy 
dependence is given by (cf. p5| for a more detailed discussion) 



ei/ \ 1 (E + 3Q) 2 
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while a transition p— >s has the energy dependence 



og, (p - s) oc • (13) 

If E <C Q the conventional energy dependence is recovered. From the above equations we 
obtain contributions to the reaction rate which are not constant (for s-wave capture) or 
proportional to Tg (for p-wave capture) in the case of small Q-values. 

Several transitions considered in this work have Q-values less than 1 MeV (e.g., all tran- 
sitions of the reaction 16 C(n,7) 17 C). For these cases the neutrons are very loosely bound and 
the bound state wave functions reach out very far. We therefore call the deviations from 
the conventional energy dependence halo effects. 

In the reaction 13 C(n,7) 14 C we also have a contribution from an incoming d-wave. The 
energy dependence of d-wave capture is, if E <C Q, proportional to E 3 ^ 2 . 

We now parameterize the total non-resonant reactions as a function of temperature Tg, 

N A < (TV > m = A + BT 9 — CTg° cnAnole"^" 1 . (14) 

The first term and the second term arise from the s-wave and p-wave contribution, respec- 
tively (cf. above). The halo effects and the d-wave contributions can be fitted as CT g D with 
only a small error. 

C. Interferences 

If the widths of the resonances are broad an interference term has to be added. The 
total cross section is then given by |26 



a(E) = a m {E) + <t t (E) + 2 [a m {E)a r {E)\ 1/2 cos[8 r {E)\ . (15) 



In this equation 5 T (E) is the resonance phase shift given by 



5 T (E) = arctan • (16) 
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Only the contributions with the same angular momentum of the incoming wave interfere 
in Eq. [1^. In the capture reactions considered in this paper we find only one case where 
the interference between a resonance and a direct capture mechanism should be taken into 
account. This is the 13 C(n,7) 14 C reaction where the p-wave resonance at 143 keV interferes 
with the p-wave contribution of the direct capture. In all other cases the interference can 
be neglected, because either the resonance is too narrow, the resonance energy is too high 
or the angular momentum of the incoming partial waves of the resonant and direct capture 
contribution differ. 

With this additional term in the cross section we have to add new terms to the reaction 
rate. We find that for our case of 13 C(n,7) 14 C the reaction rate resulting from numerical 
integration of Eq. [15] can be described by adding an interference term 

N A < av > int = ET 9 + FTg cm 3 mole~ 1 s~ 1 . (17) 



D. Nuclear Model Input Parameter 

For the calculation of the single particle amplitude in both the resonant as well as non- 
resonant neutron capture cross sections the spectroscopic factors have to be known. These 
can be obtained experimentally from single particle transfer reaction studies. For exam- 
ple, the spectroscopic factors necessary for calculating A(n,7)B can be extracted from the 
reaction A(d,p)B. The 7- widths can be extracted from reduced electromagnetic transition 
strengths. For unstable nuclei where only limited or even no experimental information is 
available, the spectroscopic factors and electromagnetic transition strengths can also be 
calculated using nuclear structure models like the shell model (SM). 

The most important ingredients in the potential model are the wave functions for the 
scattering and bound states in the entrance and exit channels. This is the case for the DC 
cross sections af c in Eq. |9| as well as for the calculation of the single-particle width r s . p . in 
Eq. [|. For the calculation of these wave functions we use real folding potentials which are 
given by |pi"l , |27 



S 



V(R) = X V F (R) = X J J p a (r 1 )p A (r 2 )i; cfr (E, p a , p A , s) dr x dr 2 , (18) 

with A being a potential strength parameter close to unity, and s — |R + r 2 — rjj, where R 
is the separation of the centers of mass of the projectile and the target nucleus. The density 
can been derived from measured charge distributions [B^] or from nuclear structure models 
(e.g. Hartree-Fock calculations) and the effective nucleon-nucleon interaction v e g has been 
taken in the DDM3Y parameterization [2~7|. The imaginary part of the potential is very 
small because of the small flux into other reaction channels and can be neglected in most 
cases involving neutron capture by neutron-rich target nuclei. 

For the calculation of the bound state wave function the parameter A is determined from 
the binding energy. In the scattering channel we try to fit A to reproduce the thermal elastic 



scattering cross section |29| . In general, we assume that A is independent of both the parity 
and the channel spin. In a few cases, however, the incoherent scattering cross section is not 
negligible. Then we need to distinguish between the different possible channel spins and 
determine A for each channel spin. One example for this is the neutron scattering on 13 C 
where different strength parameters are obtained due to the known incoherent scattering 
data. 

For the calculation of the neutron widths of resonance states the parameter A is obtained 
from the resonance energies. In a few cases the neutron widths calculated with Eq. |] can 
be compared with experimental data. This is done in Table |. In all cases the widths agree 
within a factor of 2. 

Detailed shell model calculations had to be performed to calculate the unknown excita- 
tion energies, spectroscopic factors and electromagnetic transition rates. The code OXBASH 
was used for this purpose. In most cases we need to consider a combined p- and sd-shell 
for the neutron-rich C-, N- and O-isotopes. The wave functions of these isotopes are calcu- 



lated with the interaction WBN from Ref. [31]. We included no more than lplh-excitations 



from the p- into the sd-shell. This means that for an isotope with neutron number N the 
occupation number of the sd-shell is (N — 8) for normal parity states and (N — 7) for nonnor- 
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mal parity states. For higher excitations the shell model dimensions become prohibitively 
high. 



III. NEUTRON CAPTURE BY C-, N- AND O-ISOTOPES 



Most of the neutron capture rates discussed here have been calculated previously in 
Ref. ||. However, the predicted rates were handicapped by several shortcomings, mainly in 
the determination of the nuclear structure input parameters. First, only limited experimental 
and theoretical information was available about the single particle spectroscopic factors. For 
many unstable nuclei only rather crude estimates of the spectroscopic factors were used. 
Second, the electromagnetic transition strengths were not known explicitly and systematic 
estimates were employed. Third, for some of the unstable isotopes no experimental levels 
were known above the threshold. Therefore, no resonant contributions were included for 
these reactions. 

Using the shell model we are able to calculate spectroscopic factors, electromagnetic 
transition strengths, and the resonance parameters for all reactions. While it has to be 
admitted that the reliability of the shell model calculations decreases when we approach the 
dripline, we nevertheless believe that our calculations represent an important improvement 
compared to the previous attempts. 

In the following we will discuss the neutron capture rates separately. The parameters for 
the direct capture are listed in Tables |H], UJ and |Iy]. The parameterizations of the direct 



capture contribution to the reaction rate is given in Table |V|. The resonance parameters are 



listed in Tables VI, VII and VIII. 



In Figs, m, |] and ^ we show the cross sections determined with the help of our hybrid 
model for 13 C(n,7) 14 C, 15 N(n,7) 16 N, and 18 0(n,7) 19 together with the experimental data 



from p2|-|35|j. We also compare our reaction rates with the previously proposed ones in 
Figs. |], |5] and |6|. In most cases the previous theoretical rates are taken from Ref. , except 
for the reactions 13 C(n,7) 14 C 15 N(n,7) 16 N, and 18 0(n,7) 19 where we show the comparison 
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with the experimental rates from Refs. [|32} -|3~5|1 . In Figs. |7|, [8] and ^| we compare our reaction 
rates with the reaction rates determined by the Hauser-Feshbach model [13 . 



A. 13 C(n, 7 ) 14 C 

The low energy reaction cross section of 13 C(n,7) 14 C is determined by a 2 + p-wave res- 
onance at E£ m = 143 keV which decays predominantly to the fifth excited state in 14 C at 
i? x =7. 01 Me V (0 _ ) and by the non-resonant s-wave direct capture to the ground state and 
the second excited state at E x =6.59MeV (0 + ). Additional p-wave direct capture contribu- 
tion yields from the transition to the first excited state at E x =6.09MeV (1~). Both the 
resonant cross section as well as the non-resonant cross sections have been measured re- 
cently |32],[33|. The experimental data indicate that the total cross section is dominated by 



the p-wave resonant contribution at energies above 20 keV. The non-resonant contribution 
is considerably lower and agrees well with the value extrapolated by the 1/v-law from the 
thermal cross section. 

To interpret the recent observational results in the framework of the model described 
above, we determine two potential strengths from the coherent and incoherent thermal 
scattering cross section. With these strengths we calculate the s-wave capture cross section 
in the channel spin formalism. The resulting theoretical thermal capture cross of 6.33 mb 
is larger than the experimental absorption cross section of 1.37 mb. This might partially be 
due to the El polarization induced by the incident neutron in the target nucleus |36| . In view 
of the short-range character of the nuclear forces this polarization effect is only important if 
the capture reaction takes place inside the nucleus. This is the case for the s— >p transition 
in 13 C(n,7) 14 C. Therefore we have extrapolated the experimental thermal cross section with 
an 1/f-behavior. Other direct transitions in this reaction are incoming p- waves where the 
main contributions to direct capture come from the nuclear exterior so that the polarization 
effect is supposed to be small. 

At temperatures above 0.3 GK the reaction rate is dominated by the 143 keV resonance, 
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at lower temperatures p-wave contributions determine the rate at ~ 0.1 GK and s-wave 
contributions at temperatures < 0.05 GK. Since the resonance is relatively broad we cannot 
neglect the interference between the resonance and the p-wave direct capture. The interfer- 
ence is constructive at energies lower than the resonance energy and destructive at higher 
energies. 

In Fig. the total cross section is shown with the various contributions. At higher 
temperatures the d-wave become important. The agreement with the experimental data 
from f32j|^| is satisfactory although the energy dependence of the cross section is somewhat 



different. The comparison with the rate calculated ceon the basis of experimental data 



32| , |33H is shown in Fig. [|. Experimental cross section data are only available in the energy 
range up to about 60keV. It would be interesting to measure the cross sections at higher 
energies to study the interference effects with the resonance at 143 keV and compare them 
with our calculation. The rate is considerably lower than predicted by the Hauser-Feshbach 
code SMOKER as can be seen in Fig. |7|. This can be explained by the low level density in 
the compound nucleus 14 C. 

B. 14 C(n, 7 ) 15 C 

Due to the absence of low-lying resonances the reaction rate is given by the DC con- 
tribution only. The DC transition into the ground state of 15 C is dominant. The high 
spectroscopic factor for this transition results in a clearly higher reaction rate. The reaction 
rate was calculated recently by Mengoni et al. J37|. He compares the Maxwellian averaged 
cross section with the measurement from Ref. [38] and another calculation from Ref. [[24[ . 



In Table |IX| we compare our results with these previous cross sections. Compared to 
both other calculations our result of 10.14/zb is slightly larger. All calculations are clearly 
larger than the experimental data. The reason for this discrepancy is still unknown. There 
are plans to remeasure the cross section. 
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C. 15 C(n, 7 ) 16 C 



This reaction is also purely direct since there exists no resonance level near the thresh- 



old |39[. The direct capture proceeds to the + ground state, the first excited 2 + state at 
1.766 MeV and the second excited + at 3.027 MeV. Since our spectroscopic factors - espe- 
cially for the transition to the second excited state - are higher than estimated in Ref. ||, 
the new reaction rate is larger. The new rate is also considerably lower than the predicted 
value by SMOKER (Fig. 0). 



D. 16 C(n, 7 ) 17 C 

The reaction rate is dominated by the direct transition to the second excited l/2 + state 
at 0.295 MeV. Again the spectroscopic factor and the reaction rate are higher than estimated 
by Rauscher et al. || (Fig. |j) but seems to be in good agreement with the Hauser-Feshbach 
predictions (see Fig. |?|). The contribution of the resonance at 440 keV is negligible. This 
resonance was recently identified experimentally by Raimann et al. [5tfl. The Maxwellian 



averaged cross section at 30keV was calculated by Mengoni. In Table |IX| we compare our 
results. The cross sections differ only by about 10 per cent. 



E. 17 C(n, 7 ) 18 C 

From the shell model calculations we obtain four resonances in the energy range between 
E T = 684 keV and E r = 796 keV. The reaction rate is only slighty dependent on these 
resonances, because the resonance energies are more than 700 keV above threshold For lower 
temperatures the direct capture dominates. This rate has not been calculated previously, yet 
the calculated rate is in reasonable agreement with the prediction by the Hauser-Feshbach 
model. This agreement seems to be a fortuitous circumstance. 
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F. 15 N(n, 7 ) 16 N 



The reaction has been investigated at energies between 25keV and 400 keV recently [34 
The cross section is clearly dominated by p-wave transitions for energies below 500 keV. The 
s-wave contribution, which is obtained from an extrapolation of the thermal absorption cross 
section, is negligible at thermonuclear energies. 



In Fig. ^| the cross section is compared with the measured data points from Ref. [33 
The agreement is very good. The resulting reaction rate is in good agreement with previous 
calculations ||. The present calculation is also very similar to the calculations in Ref. [33 
The small difference is due to the use of a folding potential in this work. 



G. 16 N(n, 7 ) 17 N 

Several unbound levels in 17 N are known from transfer and /3-delayed neutron decay 



studies |41]j4^] . However, no spin assignment is available for the states in the stellar energy 
range. While the shell model predicts several levels in this energy region the available 
information does not allow identification of the levels with their experimental counterparts. 
Therefore we used the shell model energies in our calculation for the resonant contribution. 
In Ref. || the experimental energy of the 198 keV resonance was used with a hypothetical 
5/2 + assignment. 

The difference in resonance energies in the previous and the present calculations explains 
the ratio of the rates. Moreover, the inclusion of the direct transition to a 5/2~ state at 
4.415 MeV with a spectroscopic factor close to unity leads to an enhancement of the reaction 
rate (see Fig. ^). Yet the rate is considerably lower than the rate predicted by the Hauser- 
Feshbach model as shown in Fig. 0. 
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H. 17 N(n, 7 ) 18 N 

Several bound states in 18 N have been identified in transfer, charge exchange, and 18 C 
/3-decay studies The measurement of the /3-delayed neutron decay of 18 C |4^| yields 



information about some neutron unbound states, but the identified levels are too high in 
excitation energy to be of relevance for the neutron capture process discussed here. Due 
to the lack of experimental information about potential resonance levels its contribution 
was calculated on the basis of shell model predictions for level energies and single particle 
strengths. The strongest contribution to its direct capture components is the transition to 
a 1~ state at 1.165 MeV with a spectroscopic factor of 0.7. This transition was not included 
in the previous calculation |§. 

Our rate is more than one order of magnitude larger than the previous estimate || but 
seems to agree reasonably well with the Hauser-Feshbach calculation. 



I. 18 N(n, 7 ) 19 N 

Only very limited experimental information is available about the level structure of 
19 N E3fl. Multi-particle transfer studies identified some of the excited bound states and 
helped to determine the corresponding excitation energies. In the study of the /3-delayed 
neutron emission of 19 C [H] several neutron unbound levels above 6.3 MeV were observed. 
No states, however, were identified near the neutron threshold of 5.32 MeV. We therefore 
rely in our estimate on shell model predictions only. A resonance a few keV above the 
threshold dominates the reaction rate for temperatures below Tg = 1. According to the 
factor exp (— 11.605-Ej/Tg) in Eq. [7] the resonant reaction rate varies by about a factor of 
2 when assuming an uncertainty of the resonance energy from 0-50 keV. According to the 
factor exp (— ll.605.Ej/Tg) in Eq. [7] the resonant reaction rate varies by about a factor of 2 
at T 9 = 1 when assuming an uncertainty of the resonance energy from 0-50 keV. 

In this temperature range our rate is therefore two to four orders of magnitude larger 
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compared to Ref. 0. For higher temperatures other resonances become important which 
is indicated by the agreement with the Hauser-Feshbach estimate for higher temperatures 
(Fig. §). 



J. 18 0(n, 7 ) 19 

The cross section for this reaction was recently measured over a wide energy range 



between 25 keV and 400 keV [ 3l| . The theoretical analysis took into account possible s- and 



p-wave transitions to the ground state and the first excited states as well as contributions 



from previously observed but unpublished higher energy resonances [45]. While the present 
evaluation is based on the same resonance parameters, the direct capture cross section has 
been reevaluated to improve the extrapolation of the data over the entire energy range. The 
present reaction rate differs only slightly from the previous result. The s-wave is determined 
by a transition to a state a few keV below the threshold. Neither the spin assignment nor 
the energy are known exactly. We extrapolate the experimental thermal cross section with 
an 1/ f-behavior to determine the s-wave contribution. In Fig. Rl we compare the calculated 



cross section with the experimental data from Ref. ||35|| . The data points at weighted neutron 



energies of 150 keV and 370 keV are clearly enhanced. 



The calculated reaction rate agrees very well with the calculation of Ref. [|35| . Again the 
small difference is primarily due to the use of a folding potential. 

K. 19 O(n, 7 ) 20 O 

Several excited states - two of them forming a doublet at 7.622 MeV excitation energy - 
are known in the compound nucleus 20 O above the neutron threshold at 7.608 MeV |f48"|,f49 



These states could be identified with shell model levels. From the shell model we also 
obtain a number of additional resonance states. The previous reaction rate included only 
two resonances. For high temperatures our rate is therefore about one order of magnitude 
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larger than the previous estimate (Fig. but seems to be in good agreement with the 
Hauser-Feshbach calculation (Fig. |]). 

L. 20 O(n, 7 ) 21 O 

The shell model predicts for the first excited state above the 5/2 + ground state of 21 
a spin and parity of l/2 + and an excitation energy of 1.33 MeV. With a high spectroscopic 
factor of 0.811 the direct p — > s transition to this state dominates the reaction rate. Pre- 
viously only the transition to the ground state was taken into account. Therefore the new 
rate is higher by a factor of approximately 5. The contributions of the two resonances are 
small. 

M. 21 0(n, 7 ) 22 

The new reaction rate is between two and three orders of magnitude higher. The spec- 
troscopic factor for the transition to the ground state is 5.222, much higher than 0.2, which 
was estimated in Ref. ||. Still, the transitions to two excited states are larger than the 
ground state transition since they are p — > s transitions. Moreover, the resonance at 14 keV 
contributes at low temperatures. 

IV. SUMMARY AND DISCUSSION 

CNO-isotopes can act as neutron poisons in s-process environments like red giants or 
asymptotic giant branch stars. They might also serve as pathways to the production of 
heavier nuclei in inhomogeneous big bang environments. In either case, high accuracy cross 
sections and reaction rates are strongly needed for relevant and precise nucleosynthesis 
predictions. 

In the present paper we combine the latest experimental information based on cross sec- 
tion measurements (for stable target nuclei) and on indirect measurements of transfer, charge 
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exchange and decay processes (for radioactive target nuclei) with theoretical direct capture 
and shell model calculations to derive reliable rates for neutron capture on neutron-rich C-, 
N-, and O-isotopes. Thus, the reaction rates provided here contain the latest information 



available and should replace the corresponding rates given in previous compilations 0,10,12 
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TABLES 

TABLE I. Comparison of the calculated neutron widths of resonance states with experimental 
values. 



reaction 


E x 




£ res [MeV] 


L n (calc) [keV] 


L n (exp) [keV] 


13 C(n, 7 ) 14 C 


8.320 


2+ 


0.143 


3.15 


3.4 ±0.6 a 


15 N(n, 7 ) 16 N 


3.360 


1+ 


0.869 


21.85 


15 ±5 b 


18 0(n, 7 ) 19 


4.109 


3/2+ 


0.152 


0.05 


< 15 c 




4.328 


5/2~ 


0.371 


6 xl0~ 4 


< 15 c 




4.582 


3/2" 


0.625 


28.5 


52 ±3 c 




4.703 


5/2+ 


0.746 


0.04 


< 15 c 



a from Ref. @ 
b from Ref. |7| 



Trom Ref. [35| 
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TABLE II. Considered transitions for the direct capture reactions on C-isotopes. Transitions 
with very small contributions are not included in the table. The spectroscopic factors are from shell 
model calculations unless stated otherwise. The s-wave transitions for the reaction 13 C(n,7) 14 C 
were obtained by extrapolating the thermal absorption cross section. Therefore the spectroscopic 
factors for these transitions are not listed. 



reaction 


Q-value (MeVJ 


J 


(MeV) 


transition 


C b 


1,j C(n,7) i4 C 


8.176 


r\-4- 

+ 


0.000 


s- 


*J-Pl/2 












d— 


*lPl/2 


1.7o4 






i — 
i 


o.uy4 


P- 


^2si /2 


U. ( OU 








6.589 


s— 


^2pi/ 2 








■5 


6.728 


P- 




0.650 






2^ 


7.012 


s— 


^ 2 P3/2 








2 


7.341 


p- 


^ld 5 /2 


0.720 


14 P|'ti n/\ 15 P 


1 QIC 
l.Zlo 


1 /9+ 
1/2 


n nnn 
u.uuu 


p- 


^2s 1/2 


u.you 






5/2+ 


0.740 


p- 


^ld 5 / 2 


0.943 


15 C(n, 7 ) 16 C 


4.251 


0+ 


0.000 


p- 


^2si /2 


0.601 






2+ 


1.766 


p- 


^ld 5 / 2 


0.493 






0+ 


3.027 


p- 


^2si/ 2 


1.344 


16 C(n, 7 ) 17 C 


0.729 


3/2+ 


0.000 


p- 


^ld 3 / 2 


0.035 






5/2+ 


0.032 


p- 


^ld 5 / 2 


0.701 






1/2+ 


0.295 


p- 


^2s 1/2 


0.644 


17 C(n, 7 ) 18 C 


4.180 


0+ 


0.000 


p- 


^ld 3 /2 


0.103 






2+ 


2.114 


p- 


^ld 5 / 2 


1.081 










p- 


^2si /2 


0.015 






2+ 


3.639 


p- 


^2s 1/2 


0.525 



from Ref. gfj 
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TABLE III. Considered transitions for the direct capture reactions on N-isotopes. Transitions 
with very small contributions are not included in the table. The spectroscopic factors are from shell 
model calculations unless stated otherwise. The s-wave transitions for the reaction 15 N(n,7) 16 N 
were obtained by extrapolating the thermal absorption cross section. These transitions are not 
listed in the table. 



reaction 


y-value (ivlev J 


J 


ti x (Mevj 


transition 




IN (^,7 ) IN 


O A Ol 

z.4»i 


z 


U.UUU 


P- 


^ld 5 /2 


u.oou 









n i on 

U.lzl) 


P" 


^2si /2 


0.460 






Q — 






P- 


1 A 

+W5/2 


n k A n a 
0.540 






i — 
1 


U.oy r 


P- 


^2si /2 


n con a 


In (n,7) in 


0.800 


5/2 


1.907 


P- 


1 A 

•>ld 5 /2 


0.207 






7 /I— 
f/Z 


o. izy 


P- 


1 A 

^ld 5 / 2 


1 /I c;7 
1.40 1 






5/2 


/I A 1 £ 

4.415 


P- 


^2si /2 


0.921 


17t\t/ \18at 

'N(n,7j °N 


2.825 


2~ 


0.121 


P- 


^ld 5 / 2 


0.700 






Q — 
O 


n 7/17 
U. 747 


P- 


^ld 5 /2 


0.689 






1 


1.165 


P- 




^2s 1/2 


0.705 


IStmY-k, -A 19 XT 

IN(n,7j In 


o.ozo 


6/ Z 


l.OOZ 


P- 


4l d 5 /2 


n £7n 
U.O f 9 






1 /9+ 
1/2 


z. 110 


s- 


^2p 3/2 


n nm 

U.UUl 










s— 


^2pi/ 2 


0.0007 






5/2- 


2.173 


P- 


^ld 5 /2 


0.354 






5/2+ 


2.375 


s— 


^2p 3 / 2 


0.001 






3/2- 


3.591 


p- 


^2si /2 


0.461 






3/2+ 


3.799 


s— 


^2p 3/2 


0.007 










s— 


^2pi/ 2 


0.001 






1/2- 


4.126 


p- 


^2s 1/2 


0.606 






3/2- 


4.438 


p- 


^2si /2 


0.035 






1/2- 


5.101 


p- 


^2si /2 


0.440 






5/2+ 


5.130 


s— 


+2p 3 /2 


0.001 






3/2- 


5.173 


p- 


^2s 1/2 


0.097 



from Ref. (5(J 
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TABLE IV. Considered transitions for the direct capture reactions on O-isotopes. Transitions 
with very small contributions are not included in the table. The spectroscopic factors are from shell 
model calculations unless stated otherwise. The s-wave transitions for the reaction 18 0(n,7) 19 
were obtained by extrapolating the thermal absorption cross section. Therefore the spectroscopic 
factors for these transitions are not listed. 



reaction 


y-value (IVlev J 


J 


h x (JVleV) 


transition 


/~>2 c 
G D 


i0 0(n,7j iy O 


3.957 


5/2 + 


0.000 


P- 




0.570 a 






l/2 + 


1.472 


P- 


-^-iSi/2 


1.000 a 






/9 — 
Of A 


"i Q/l fi 


s— 


^2p 3/2 




19 O(n, 7 ) 20 O 


7.606 


0+ 


0.000 


P- 




3.427 






2+ 


1.674 


P- 


^ld 5 / 2 


0.731 










P- 


^2s 1/2 


0.142 






4+ 


3.570 


P- 


^ld 5 / 2 


1.021 






2+ 


4.072 


P- 


^2s 1/2 


0.573 






3+ 


5.447 


P- 


^2s 1/2 


0.817 


20 O(n, 7 ) 21 O 


3.806 


5/2+ 


0.000 


P- 


^ld 5 / 2 


0.345 






1/2+ 


1.330 


P- 


^2s 1/2 


0.811 


21 0(n, 7 ) 22 


6.850 


0+ 


0.000 


P- 


^ld 5 / 2 


5.222 






2+ 


3.374 


P- 


^2s 1/2 


0.822 






3+ 


4.830 


P- 


^2s 1/2 


0.771 



a from Ref. @ 
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TABLE V. Parameters for the direct capture reactions. The interference term with the param- 
eters E and F is only used for the reaction 13 C(n,7) 14 C. 





A 


B 


C 


D 


E 


F 


13 C(n, 7 ) 14 C 


182.310 


3296.558 


-5534.721 


2.009 


1794.927 


-1103.269 


14 C(n, 7 ) 15 C 


— 


4754.286 


752.370 


1.630 






15 C(n, 7 ) 16 C 


— 


2637.750 


304.878 


1.644 






16 C(n, 7 ) 17 C 


— 


2861.539 


1166.516 


1.311 






17 C(n, 7 ) 18 C 


— 


1334.578 


337.957 


1.472 






15 N(n, 7 ) 16 N 


3.18 


3783.415 


335.198 


1.716 






16 N(n, 7 ) 17 N 


— 


3649.913 


437.549 


1.633 






17 N(n, 7 ) 18 N 


— 


3417.690 


358.029 


1.660 






18 N(n, 7 ) 19 N 


13.838 


4051.118 


966.727 


1.412 






18 0(n, 7 ) 19 


21.357 


8300.018 


596.897 


1.770 






19 O(n, 7 ) 20 O 




7275.368 


432.266 


1.747 






20 O(n, 7 ) 21 O 




6474.727 


493.516 


1.750 






21 0(n, 7 ) 22 




7327.938 


543.151 


1.747 







23 



TABLE VI. Adopted values for the resonance parameters for capture reactions on carbon 
isotopes. The neutron and 7 widths are calculated as explained in the text unless stated otherwise. 



Reaction 


E x 




Exes 


r n 


r 7 






[MeV] 




[MeV] 


[eV] 


[eV] 


[eV] 


13 C(n, 7 ) 14 C 


8.320 


2+ 


0.143 


3400 a 


0.215 b 


0.269 


16 C(n, 7 ) 17 C 


1.180 


1/2" 


0.451 


950 


2.82 -10~ 4 


2.82 -10~ 4 


17 C(n, 7 ) 18 C 


4.864 


4+ 


0.684 


50 


6.11 -10~ 4 


6.87 -10" 4 




4.915 


3+ 


0.735 


8 450 


7.37 -10~ 3 


6.45 -10~ 3 




4.972 


l - 


0.792 


5 440 


5.85 -10~ 2 


0.022 




4.976 


2+ 


0.796 


8 510 


1.85 -10~ 3 


1.16 -10~ 3 



a from Ref. @ 
b from Ref. ^2| 



24 



TABLE VII. Adopted values for the resonance parameters for capture reactions on nitrogen 
isotopes. The neutron and 7 widths are calculated as explained in the text unless stated otherwise. 



Reaction 


E x 


,r 


E 


r„ 


r 7 


W7 




[McV] 




[McV] 


[eV] 


[eV] 


[eV] 


15 N(n, 7 ) 16 N 


3.360 


1+ 


0.862 


15 000 a 


0.455 


0.341 


16 N(n, 7 ) 17 N 


5.904 


7/2- 


0.021 


0.032 


4.80 -10~ 2 


0.015 




6.121 


5/2+ 


0.238 


1.2 


4.80 -10~ 2 


0.027 




6.325 


3/2+ 


0.442 


20 


5.46 -10~ 2 


0.022 




6.372 


7/2+ 


0.489 


20 


1.52 -10 -2 


0.012 




6.373 


5/2+ 


0.490 


600 


0.110 


0.066 




6.470 


1/2+ 


0.587 


1750 


2.510 


0.501 




6.685 


3/2- 


0.802 


12 500 


5.660 


2.263 




6.737 


7/2+ 


0.854 


70 


4.17 -10- 2 


0.033 




6.835 


3/2+ 


0.952 


360 


0.478 


0.191 


17 N(n, 7 ) 18 N 


2.875 


3- 


0.050 


1.6 


9.61 -10~ 3 


0.017 




2.949 


2+ 


0.124 


1390 


9.72 -10~ 2 


0.121 




3.068 


1+ 


0.243 


1060 


0.209 


0.157 




3.374 


3+ 


0.549 


50 


0.107 


0.187 




3.437 


V 


0.612 


10180 


0.517 


0.646 




3.631 


0+ 


0.806 


2 350 


6.26 -10~ 2 


0.016 




3.644 


l - 


0.819 


2 370 


0.319 


0.239 




3.722 


2+ 


0.897 


190 


4.48 -10~ 2 


0.056 


18 N(n, 7 ) 19 N 


5.335 


3/2+ 


0.007 


7.8 


0.128 


0.050 




5.479 


3/2+ 


0.151 


70 


0.072 


0.029 




5.498 


7/2+ 


0.170 


40 


2.40 -10~ 2 


0.019 




5.634 


9/2" 


0.306 


120 


5.72 -10~ 4 


5.7 -10~ 4 




5.770 


5/2+ 


0.442 


170 


4.28 -10~ 2 


0.026 




5.778 


3/2+ 


0.450 


120 


0.320 


0.128 




5.858 


3/2- 


0.530 


8 470 


2.520 


1.008 




5.955 


7/2+ 


0.627 


180 


8.26 -10- 2 


0.066 




6.006 


5/2+ 


0.678 


190 


4.85 -10- 2 


0.029 




6.094 


7/2- 


0.766 


9 990 


0.306 


0.245 




6.125 


7/2+ 


0.797 


950 


5.08 -10~ 2 


0.041 




6.128 


3/2+ 


0.800 


140 


1.130 


0.448 




6.130 


5/2+ 


0.802 


310 


7.86 -10~ 2 


0.047 




6.152 


1/2+ 


0.824 


2 090 


1.100 


0.220 
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TABLE VIII. Adopted values for the resonance parameters for capture reactions on oxygen 
isotopes. The neutron and 7 widths are calculated as explained in the text unless stated otherwise. 



Reaction 




r 


-E'res 


r n 


r 7 


w 7 






[MeV] 




[MeV] 


[eV] 


[eV] 


[eV] 




18 0(n, 7 ) 19 


4.109 


3/2+ 


0.152 


50 


1.4 -10" 2 


0.028 






4.328 


5/2" 


0.371 


0.6 


0.600 


0.900 






4.582 


3/2" 


0.625 


52 000 a 


1.900 


3.800 






4.703 


5/2+ 


0.746 


40 


0.500 


1.483 




19 O(n, 7 ) 20 O 


7.622 


3- 


0.015 


11 


6.52 -10~ 2 


0.038 






7.622 


4+ 


0.015 


1.1 


1.11 • 10~ 2 


8.24-10" 


-3 




7.638 


4- 


0.031 


2.3 


1.82 -10~ 2 


0.014 






7.646 


2~ 


0.039 


3.3 


0.109 


0.044 






7.739 


1- 


0.132 


550 


12.33 


3.015 






7.754 


4+ 


0.147 


1190 


0.301 


0.226 






7.855 


5- 


0.248 


1.8 


1.45 -10~ 2 


0.013 






7.970 


2+ 


0.363 


2 480 


4.46 -10~ 2 


0.019 






8.160 


3- 


0.553 


2 890 


0.352 


0.205 






8.403 


1- 


0.796 


9150 


0.362 


0.090 






8.439 


2 - 


0.832 


2090 


0.333 


0.139 






8.533 


3" 


0.926 


8 840 


0.688 


0.401 






8.552 


2 - 


0.945 


1480 


0.238 


0.099 






8.558 


3+ 


0.951 


20 000 


0.320 


0.187 




20 O(n, 7 ) 21 O 


4.343 


5/2" 


0.536 


0.8 


0.121 


0.316 






4.765 


3/2 - 


0.958 


5 450 


0.631 


1.262 




21 0(n, 7 ) 22 


6.863 


4+ 


0.014 


0.7 


6.83 -10~ 4 


5.12-10" 


-4 




7.357 


4+ 


0.508 


26 810 


2.20 -lO" 3 


1.65-10" 


-3 




7.397 


2 - 


0.548 


780 


0.014 


5.83-10" 


-3 




7.472 


3" 


0.623 


790 


0.023 


0.013 






7.515 


1- 


0.666 


42 710 


0.146 


0.036 





'from Ref. || 
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TABLE IX. Comparison of the Maxwellian averaged cross section < av > jkT at kT = 30keV 
with other works. 



reaction 
14 C(n, 7 ) 15 C 
16 C(n, 7 ) 17 C 



this work 
10.14 
4.71 



Ref. |7| 
8.3 
4.3 



Ref. H 
8.4 



Ref. H 
1.87 ± 0.43 
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FIGURES 

FIG. 1. Calculated cross section of the reaction 13 C (11,7) 14 C compared with experimental data 
from Ref. pLB3]. 



FIG. 2. Calculated cross section of the reaction 15 N(n,7) 16 N compared with experimental data 



from Ref. [34|. 



FIG. 3. Calculated cross section of the reaction 18 0(n,7) 19 compared with experimental data 



from Ref. [35 



FIG. 4. Comparison of our new reaction rates for neutron capture on C-isotopes with previous 
direct capture calculations. Shown is the ratio of the new rates to the rates published in Ref. |J 
for all rates except 13 C(n,7) 14 C where the ratio to the rate determined from the experimental data 
of Ref. HI! is shown. 



FIG. 5. Comparison of our new reaction rates for neutron capture on N-isotopes with previous 
direct capture calculations. Shown is the ratio of the new rates to the rates published in Ref. |J 
for all rates except 15 N(n,7) 16 N where the ratio to the rate determined from the experimental data 



of Ref. [34] is shown. 



FIG. 6. Comparison of our new reaction rates for neutron capture on O-isotopes with previous 
direct capture calculations. Shown is the ratio of the new rates to the rates published in Ref. || 
for all rates except 18 0(n,7) 19 where the ratio with the rates determined from the experimental 



data of Ref. [35] is shown. 



FIG. 7. Comparison of our reaction rates for neutron capture on C-isotopes with 
Hauser-Feshbach calculations obtained with the code SMOKER. 



FIG. 8. Same as Fig. [?] for neutron capture of N-isotopes. 



FIG. 9. Same as Fig. ^ for neutron capture of O-isotopes. 



33 



p-wave + | 
resonance , 



1 00 

a(nb) 

1 



1 



15 N(n, 7 ) 16 N 
16 N(n, 7 ) 17 N 
17 N(n, 7 ) 18 N 
18 N(n, 7 ) 19 N 



18 0(n, 7 ) 19 
19 O(n, 7 ) 20 O 
20 O(n, 7 ) 21 O 
21 0(n, 7 ) 22 



